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Abstract 

We review recent experimental measurements of the Fermi surface of the iron-pnictide superconductor 
LaFePO using quantum oscillation techniques. These studies show that the Fermi surface topology is 
close to that predicted by first principles density functional theory calculations, consisting of quasi-two- 
dimensional electron-like and hole-like sheets. The total volume of the two hole sheets is almost equal 
to that of the two electron sheets, and the hole and electron Fermi surface sheets are close to a nesting 
condition. No evidence for the predicted three dimensional pocket arising from the Fe band is found. 
Measurements of the effective mass suggest a renormalisation of around two, close to the value for the overall 
band renormalisation found in recent angle resolved photoemission measurements. 
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1. Introduction 

Experimental determinations of the Fermi sur- 
face are important for a number of reasons. Firstly, 
although calculations of the band structure using 
density functional theory within the local density 
approximation (and its variants) give a good gen- 
eral picture, in many cases the calculated band po- 
sitions are found to be accurate to only a few tens 
of meV [1, 2, 3]. These relatively small errors can 
lead to substantial differences in the size and shape 
of the Fermi surface. For the iron-pnictides this will 
have a bearing on, for example, how close to perfect 
nesting the various Fermi sheets are and also the 
magnitude and anisotropy of the Fermi velocities. 
This will in turn lead to corrections to calculations 
of the magnetic excitation spectrum and the inter- 
actions that may be responsible for superconduc- 
tivity, plus other physical properties, such as Hc2 
and its anisotropy. Secondly, conventional band 
structure calculations only include many-body ef- 
fects at the mean field level, and so in strongly inter- 
acting systems there can be substantial renormal- 
isation of the effective masses or Fermi velocities 
from electron-phonon or electron-electron interac- 
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tions. These mass renormalisations need to be de- 
termined experimentally. Thirdly, in systems with 
structural or magnetic instabilities (such as charge 
density waves and spin density waves) there may be 
subtle changes in the lattice structure which could 
lead to a reconstruction of the Fermi surface, pos- 
sibly breaking up a large Fermi surface into many 
small pockets. 

Magnetic field induced quantum oscillations ef- 
fects in, for example, the magnetisation [the de 
Haas-van Alphen effect (dHvA)] or the resistance 
[Shubnikov-de Haas effect (SdH)] have been used 
for more than half a century to determine the 
Fermi surfaces of metals [3]. The power of the 
technique is that it provides very accurate mea- 
surements of the Fermi surface cross-sectional ar- 
eas (typically with an accuracy better than 0.1% 
of the area of the Brillouin zone), effective masses 
(typically with an accuracy better than ~ 5%), and 
quasiparticle scattering rates. There are, however, 
several experimental constraints. The most seri- 
ous is that the dHvA signal is damped exponen- 
tially by impurity scattering, which broadens the 
Landau levels and therefore reduces the size of the 
anomaly which occurs when they cross the Fermi 
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level. This reduction is given by the Dingle factor, 
Rd — exp(— 7rfci?/ei?£), where kp is the orbitally 
averaged Fermi wavevector for the particular Fermi 
surface orbit, and £ is the inelastic quasiparticle 
mcan-free-path (which is often shorter than the cor- 
responding quantity derived from transport proper- 
ties as low angle scattering events are not weighted 
out). This means that dHvA oscillations are only 
observable at high magnetic field in very pure sam- 
ples. Finite temperature also reduces the size of 
the oscillatory signal due to the thermal smearing 
of the discontinuity in occupation number at the 
Fermi level. The signal is reduced by a factor 2 
when T ~ 0.12 B/m* (m* is the quasiparticle ef- 
fective mass measured in units of the free electron 
mass TOe). For an effective mass of ~ 2 and mag- 
netic field of 15 T the amplitude of the signal is 
reduced by a factor ten from its zero temperature 
value at T ~ 2 K, with a further attenuation of 
approximately one order of magnitude per Kelvin 
at higher temperature. Performing experiments at 
such low temperatures means that for most pnictide 
materials we are usually in the regime where the 
upper critical field is very high (>60T). Although 
it is possible to observe dHvA oscillations in the 
superconducting mixed state, the signal 

suffers further attenuation as the field is reduced 
below Hc2, 01 hence oscillations are typically only 
observed near to Hc2- 

A complication with understanding dHvA results 
is that it is not possible to determine experimentally 
the fc-space location of the observed Fermi surface 
orbits, and so interpretation of the data usually re- 
lies on a comparison with band-structure calcula- 
tions. Despite these limitations, quantum oscilla- 
tions are the most precise technique to determine 
Fermi surface geometry, including full information 
about the three dimensional topology. 

Another technique which has been used exten- 
sively in recent years to determine Fermi surfaces 
is angle resolved photoemission (ARPES). The 
ARPES technique has moderate resolution of the 
in-plane momentum but can provide a direct k- 
resolved map of the Fermi surface The kz mo- 
mentum cannot be probed directly although the 
kz point at which the in-plane momentum is mea- 
sured may be controlled to some extent by varying 
the photon energy. Unlike dHvA which is a bulk 
probe, ARPES probes just a few unit cells below 
the surface. The termination of the surface can 
be an issue, especially if the compound does not 
cleave on an electrically neutral layer . ARPES 



and dHvA are therefore complementary techniques, 
which both have their strengths and weaknesses. 

At the time of writing, dHvA experiments have 
been reported in just two iron-pnictide materials. 



Superconducting LaFcPO [11|, which will be 
the focus of this article and undoped SrFc2As2 fl3|- 
This latter compound is non-supcrconducting and 
has a antiferromagnetic ground state - which makes 
the Fermi surface reconstruct into many small pock- 
ets. 

The iron phosphide superconductor LaFePO, is 
isostructural with superconducting LaFeAsOi-^^, 
and the two materials are predicted to have very 
similar electronic structure. Its relatively low Tc ~ 
6 K has been linked to the fact that the Fe- 
P bond angles depart substantially from those of a 
regular tetrahedron [l^l . There is some controversy 
about whether stoichiometric LaFePO is supercon- 
ducting, or whether a few percent of oxygen va- 
cancies are present even in the nominally undoped 
material Oxygen deficiencies of around 1% 

are difficult to determine on small single crystals 
and so for the moment this remains an unresolved 
question. One significant difference from the ar- 
senide compounds is that LaFePO is not magnetic 
even in its non-superconducting state [l3| and so 
we do not expect any Fermi surface reconstruction. 
Another difference is that recent penetration depth 
measurements suggest a superconducting order 
parameter with line nodes, rather than the fully 
gapped state that has been observed in several As- 
based materials 1^ 23|- This may mean that the 
order parameter symmetry in iron-pnictides may 
vary between materials. A detailed knowledge of 
the differences in electronic structure may be im- 
portant in understanding how this arises. 

The low Tc of LaFePO leads to a correspond- 
ingly low value of ^oHc2 — 0.7 T, for H\\c. This, 
combined with the very high purity levels found for 
flux grown crystals (residual resistance ratios up to 
~ 90) makes LaFePO an excellent choice for dHvA 
studies. To date, two groups have reported dHvA 
results for LaFePO. This article will focus mainly 
on our own data ,111] but will include comparisons 
to the data of Sugawara et al. 1'^ where appropri- 
ate. 

2. Experimental details 

2.1. Sample growth and characterization 

Single crystals of LaFePO were grown from a 
molten Sn flux. La, Fe203, P and Sn mixed in the 
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Figure 1: The temperature dependence of the interlay er re- 
sistivity of LaFePO (sample C). 



molar ratio 3:1:2:24, were placed in alumina cru- 
cibles and sealed in quartz tubes under a small par- 
tial pressure of argon. Typically the melt was held 
at 1190° C for IShrs, before cooling at 10°C/hour 
to 650°C. The remaining liquid was decanted at 
the base temperature using a centrifuge, and the 
crystals removed from the crucible. Excess Sn on 
the surface of the crystals was removed by etch- 
ing in dilute HCl followed by an immediate rinse 
in methanol. Measurements of the magnetic pen- 
etration depth [is'l showed no evidence for any Sn 
flux remaining on the samples (with a resolution 
of about 1 iiir?). Crystals of up to 0.4 mm on a 
side were produced by this process although, as all 
our dHvA measurements were performed using the 
micro-cantilever method, only small single crystals 
were selected. It is likely that the smaller crys- 
tals would be more homogeneous than the larger 
ones. Some growth runs were made by an alterna- 
tive route using La203 instead of Fe203 to intro- 
duce oxygen into the melt. Details are described in 
Ref. [16|. Sugawara et al. used a very similar 
Fe2 03 based sample preparation route. 

The crystals were characterized by x-ray diffrac- 
tion, electrical resistivity and specific heat. The x- 
ray results show that crystals are tetragonal, space 
group Pnmm, with cell dimensions a = 3.941(2)A, c 
= 8.507(5)A and atomic positions .Zp=0. 63477 and 
ZLa= 0.1489 in agreement with other previous re- 
sults [l3|. The crystallographic orientation of the 
crystals was identified using Laue diffraction. 

Resistivity measurements were conducted using a 
standard low frequency ac, 4-probe technique, with 
contacts applied with Epotek H20E paste. Fig. [T] 



Figure 2: Specific heat anomaly at the superconducting tran- 
sition in a single crystalline sample of LaFePO (sample A). 
AC = C{0) — C(B), i.e., the zero field specific heat with nor- 
mal state background subtracted by using data taken in a 
field large enough to suppress superconductivity (B = IT). 



shows the c-axis resistivity of sample C, which in- 
dicates conventional metallic behaviour, supercon- 
ductivity at low temperature, and a residual resis- 
tance ratio R{mOK)/R{Tc) ~ 52. Similar resis- 
tance data are also reported by Sugawara et al. 
[l^ . In-plane resistivity has also been measured 
on other samples and the anisotropy pd Pa, which 
was relatively independent of temperature, was in 
the range 13 to 17 [11]. 

Heat capacity was measured using an ac tech- 
nique, with modulated light as a heating source 
and a 12/im thermocouple as a thermometer. This 
allowed us to measure the same small single crys- 
tals that were used for the dHvA study (typical 
dimensions 100 x 100 x 2Qpvo?). The data shown in 
Fig. m shows a sharp, well defined anomaly at the 
superconducting transition (Tc midpoint = 5.8 K) 
showing that there is bulk superconductivity in our 
LaFePO samples. 

2.2. Micro- cantilever torque measurements 

Due to the small size of the best crystals of 
LaFePO, the most suitable techniques for mea- 
suring quantum oscillations are micro-cantilever 
torque or resistivity. For the torque measurements 
the samples were attached with quick dry epoxy to 
the end of the self-sensing microcantilever (Seiko 
instruments model: SSI-SS-ML-PRC120). The re- 
sistance of the cantilever is measured with a conven- 
tional Wheatstone bridge arrangement, using an ac 
current excitation of typically / — 20/iA at a fre- 
quency of 72 Hz. The cantilever was mounted on 
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Figure 3: a) Torque measurements on LaFePO measured at 
T = 0.35 K and in magnetic fields up to 18 T for different 
magnetic field directions (sample B). The arrow indicates 
the position of Hc2 ■ b) Torque measurements up to 45 T for 
sample A. c) Field dependence of the interlayer resistivity in 
magnetic fields up to 45 T for sample C. 



a single axis rotator and placed in a '^He cryostat, 
in the bore of a 19 T superconducting magnet (in 
Bristol) or the 45 T hybrid magnet in Tallahassee, 
Florida. We have used this technique previously for 
studies of MgB2 i21] and cuprate superconductors 

Sugawara et al. |12l | used the same cantilevers, 
but with a field modulation technique; a dc current 
is passed through the cantilever and an additional 
oscillatory field {Bac — lOmT, / =lf Hz) is applied 
yielding voltage oscillations directly related to the 
torque. 



temperatures T « 0.35(5) K on different samples of 
LaFePO. In low magnetic fields the torque signal 
is characteristic of a type II superconductor in the 
vortex state with ^^0-^12 ~ 0.68 T (Fig. Hi). As the 
system is rotated away from H\\c the position of 
Hc2 shifts to much higher values {fioH^'^ « 7.2 T). 
These are the characteristics of an anisotropic su- 
perconductor with an anisotropy of 7^ « 10 [j^ ]. 
similar to that found in SmFeAsOi-fF^ 1251 but 
much larger than that for (Ba,K)Fe2As2 [26|. The 
torque response in LaFePO is mostly reversible for 
field sweeps up and down suggesting weak vortex 
pinning in our crystals. 

For magnetic fields higher than Hc2 the back- 
ground torque signal varies quadratically with H 
as expected for a simple paramagnet. Above 9 T 
the quantum oscillations start to become visible su- 
perimposed on the background signal. The inter- 
layer resistance shows a positive magnetoresistance 
and above 30 T a weak oscillatory signal is visi- 
ble (Fig. [St). By subtracting the background signal 
from the torque and resistivity data using a polyno- 
mial (5**^ order) the richness of quantum oscillation 
signal is clearly visible (Fig.[4|). 

The fast Fourier transform (FFT) of the oscil- 
latory data for different samples is shown in Fig. 
[S] We identify the dHvA frequencies F, which 
are related to the extremal cross-sectional areas 
Ak, of the FS orbits via the Onsager relation, 
F = {h/2TTe)Ak as the peak in the FFT. The an- 
gular dependence of these frequencies as the mag- 
netic field angle, 9, is rotated away the c axis 
(9 = 0°) towards the conducting planes {9 = 90°) 
[2 7| allows us to construct a three-dimensional pic- 
ture of the shape and size of the Fermi surface. 
Fig. shows such a detailed angular dependence 
for sample B, which was measured up to 18 T at 
T =0.3(1) K. In this particular experiment four dif- 
ferent oscillatory components can be seen, whose 
frequency and amplitude is strongly dependent on 
angle. The frequencies shift to higher values as 9 
is increased. This is described approximately by 
F cx 1/ cos 9, strongly suggesting that the signals 
arise from orbits on quasi-two-dimensional Fermi 
surfaces sheets. The amplitude of the oscillatory 
part of torque. 



3. dHvA results 

Fig. [31 shows raw torque and resistivity measure- 
ments in high magnetic fields up to 45 T and at low 
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(see Ref. [J]), is determined by the Fermi surface 
anisotropy and thus vanishes near symmetry axes 
and for spherically symmetric Fermi surfaces. 
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Figure 6; a) Angle dependence of all observed frequencies. Different symbols correspond to: sample A (open circles), sample B 
(filled circles), sample C (diamonds) and sample D (stars). Possible harmonics of the main frequencies are shown by triangles 
and the dashed lines indicate their calculated location. Solid lines overlapping the 7, <5 and e frequencies are fits to a 1/cos 6 
dependence, b) Angle dependence of FcosB. Solid lines are calculations for a simple cosine warped cylinder (see text). The 
position of the 'Yamaji angle' is indicated by the arrow. Note the logarithmic vertical axis on both parts of the figure. 



Table 1: Experimental effective masses (m*) and frequencies from dHvA data for samples A, B and D (our results) and sample 
S (after Ref.O)^ 
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Four frequencies have significantly larger ampli- we label ai and a2, have very similar frequencies, 
tudes than the others. The first two of these, which with values near F ^ 1 kT {AF « 35 T). Near H||c 
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Figure 4: a) Oscillatory part of the torque (dHvA) (sample 
A) and b) resistance (SdH) (sample C) in magnetic fields up 
to 45 T for different magnetic field angle, 8, with respect to 
the c axis. 



these frequencies are separated enough for distinct 
peaks to be observed directly in the FFT. At higher 
angles these peaks tend to merge, but beating can 
be seen in the raw data, indicating the survival of 
two components, whose frequencies can be found 
by fitting the torque data directly. Two higher fre- 
quencies, Pi (« 1.9 kT) and P2 {- 2.4 kT), also 
have a rather strong amplitude and have frequen- 
cies that rise rapidly as 9 is increased, crossing near 
e = 45°. 

Frequencies ai,a2,(3i and P2 were observed in 
all measured samples. A, B, D and L, as shown in 
Fig. [5] Sugawara et al. [ij], also observed these 
four orbits with almost identical frequencies and 
angle dependence. 
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Figure 5; a) The angle dependence of the Fourier transform 
spectra on LaFePO measured at T = 0.35 K for the field 
range 10-18 T (sample B). The large peak near 45° is due 
to the enhancement of the dHvA at the Yamaji angle (see 
text). Different runs are shifted for clarity and the angle 
step corresponds to Ri 3° . b) Fourier transform spectra for 
different field orientations for sample A (field range 15-45 T), 
sample L (field range 15-50 T) and sample D (field range 10- 
18 T). The positions of the orbits with smaller amplitudes, 
7, 5, and e are shown by arrows. The positions of the main 
frequencies corresponding to orbits near extremal areas of 
Fermi surface of LaFePO are shown by their labels ai, 0.2, 
/3l and /32. 



In further experiments additional frequencies 
were observed. Sample D was studied in fields up 
to 18 T, the same range as sample B, but comes 
from the batch with a higher residual resistivity ra- 
tio (~ 85). In this crystal we observe an additional 
orbit, denoted 7, with F « 0.7 kT which has a sim- 



ilar amplitude to the (32 orbit. Samples A and L 
were measured in much higher fields (up to 50 T) 
and there we observed in addition to a, /? and 7 
a further two distinct frequencies which we label 5 
and e. These were not observed by Sugawara et 
al. [l^ . probably because of the lower maximum 
field available in their experiment. The values of all 
the observed frequencies, extrapolated to H||c, are 
listed in Table[TJ We observe a very good agreement 
for all the measured samples, with the frequencies 
agreeing to within a few percent. 

The angular dependence of the frequencies as the 
magnetic field is tilted away from the c axis is shown 
in Fig. [SI For a two-dimensional cylindrical Fermi 
surface the dHvA frequency varies like 1 / cos 6 and 
deviations from this indicate the degree of warping. 
This is shown in Fig.jljD by plotting F{9) cos 9. The 
Fermi surfaces of branches ai, a2 and 7 are close to 
perfect cylinders as well as those for 5 and e (within 
the limitation of the available data). The upward 
curvature on the F cos 9 plot of the /3i data shows 
that this orbit originates from a section of Fermi 
surface that is locally concave along the c-axis (e.g., 
the minimum of a cosine warped tube), and likewise 
(32 is from a convex (maximum) section. 

For a quasi-two-dimensional Fermi surface with 
weak cosine c-axis warping, Yamaji has shown 
that the angular dependence of the two frequencies, 
F-t, arising from extremal areas of a single warped 
cylinder is given by, 

F±(9) cos 9 = {h/2Tre)[Kkp ± inmetcJoickp tan^)]. 

In this expression, the first term is the mean fre- 
quency {h/2TTe)nkl = (F0i„ + F0,,)/2, the prefac- 
tor of the second term {h/2Tre)4TTmetc is (f,°ax ~ 
-^min)/2 and Jo is the Bessel function (c is the lattice 
parameter and tc parameterises the warping). Here 
i^^ijj and -F^ax represent minimum and maximum 
frequencies at 9 = 0, and may be taken directly 
from the data, so this equation has no free parame- 
ters. A plot of this function is shown in Fig.jSb, and 
provides a very good description of the data. For 
certain angles, 9, where the Bessel function Jq = 
the two frequencies cross. At this point all the cross 
section areas are equal and their contributions to 
the total magnetization interfere constructively to 
give a sharp maximum in the amplitude of the sig- 
nal (peak effect). This 'Yamaji angle' is observed 
close to 6* - 45(2)° for the /3 orbits in LaFePO, 
(see Fig. Eb) and the peak effect is clearly seen in 
Fig. [5^. These observations show convincingly that 



(3i and P2 are the minimum and maximum frequen- 
cies of the same quasi-two-dimensional Fermi sur- 
face sheet. A similar behaviour, albeit with much 
smaller warping, is observed for the a frequencies 
indicating that these too originate from a single 
quasi-two-dimensional Fermi surface sheet. 

3.1. Torque interaction 

The total torque response, both background and 
oscillatory, causes a deflection of the lever, and so 
the angle 9 is not quite constant during the field 
sweep. This effect produces harmonics of the main 
frequencies (indicated in Fig. [5^) and may also 
produce mixing of the main frequencies, such as 
/3i_2 i ai,2, 2a, 2/3, etc. and so it is important to 
identify any such frequencies. Close to 6* = we 
find that the frequency of S (Fg) is very close to 
— Fa^ ■ However, if 5 was produced by torque 
interaction we would expect it to have an angle de- 
pendence intermediate between (32 and a2. This 
does not appear to be the case (see Fig. [Sb), how- 
ever more detailed data will be required to settle the 
issue. The 7 and e frequencies cannot be generated 
by this effect. 

3.2. Effective mass 

The cyclotron effective mass of the quasiparti- 
cles on the various orbits were determined by fitting 
the temperature dependent amplitude of the oscilla- 
tions to Rt = X/ sinh(X) where X = 14.69m*T/B 
[J] , as shown in Fig. [71 These measurements were 
made with the field close to the c-axis and the value 
for 9 = Q° was inferred using the usual approxima- 
tion that the mass scales with the dHvA frequency 
(for sample B the mass was measured at several dif- 
ferent angles to confirm this). The obtained masses 
range between 1.7 me and 2.1 Trig and are listed in 
Table [TJ These values are slightly smaller than the 
those reported in Ref. [l3|, but within the quoted 
errors. 

From the field dependence of the amplitude at 
constant temperature we can estimate the inelas- 
tic quasiparticle mean-free-path, ^, from the Din- 
gle factor Ro. We find a mean free path of ~ 
1300 A and 800 A for a and (3 orbits, respec- 
tively. A similar value of 940 A was found for the 
(3i orbit by Sugawara et al.\Y^. For the 5, 7 and e 
orbits the scattering rate is a factor ~ 2 — 3 larger 
than that for the a orbits although the small size 
of the signal precludes any detailed analysis. 
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Figure 8: The calculated band-structure of LaFePO (without spin-orbit corrections) along principle symmetry directions. The 
colour scale reflects the total Fe d character. Bands with Fe d^2 character greater than 50% are coloured black. The inset 
shows the Brillouin zone with symmetry labels. 



4. Band structure 

The band structure of LaFePO has been calcu- 
lated by Lebegue [2£|] , and consists mainly of tubu- 
lar two dimensional sheets of Fermi surface; two 
hole like tubes at the zone centre and two electron- 
like tubes at the zone corner. In addition there is a 
more three dimensional hole-like pocket near the Z 
point (see Fig. [5]). 

We have repeated these calculations so that we 
can calculate the dHvA parameters and make a de- 
tailed comparison to our data. Our own calcula- 
tions were carried out using the WIEN2K pack- 
age (30l |. which is an implementation of a full- 
potential, augmented-plane-wave plus local-orbital 
scheme. We used a generalized-gradient approxi- 
mation form for the exchange correlation potential 
[sH and the experimental crystal structure. We 
used 10^ fc-points in the full Brillouin zone for con- 
vergence and generally 10^ /c-points for the Fermi 
surface and dHvA parameter calculations. Larger 
fc-point densities were used for the two-dimensional 
cuts (typically 10^ fc-points in each 2D cut). The 
calculations were performed with and without the 



spin-orbit correction which was found to produce 
small but important differences in the band struc- 
ture. 

The band structure without spin-orbit correc- 
tions is shown in Fig. [S] These results are in good 
agreement with those of Lebegue [i^. The bands 
crossing the Fermi level are mostly of Fe d character 
as emphasized in the plot, and mostly weakly dis- 
persive along the c-direction (see the sections of the 
plot r ^ Z, and M). The exception is a band 
which has mostly Fe d^^ character which is empha- 
sized in black in the plot. In Fig. [9] we show the 
resulting Fermi surface and in Fig. [10] two dimen- 
sional slices of the Fermi surface emphasizing the 
character of the bands. In Fig.[Tn]the plot is colour 
coded according to the strongest band character at 
each individual fc-point. At the zone corner there 
are two electron-like quasi-two-dimensional Fermi 
surface sheets, which touch along the (100) direc- 
tions. In the Brillouin zone corresponding the iron 
sublattice (which is \/2 larger and rotated by 45° 
compared to the actual Brillouin zone correspond- 
ing to the real crystal structure [s^]), these sheets 
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Figure 7; The temperature dependence of the Fourier am- 
phtude for sample A {9 = 2° and field window between 25 T 
and 45 T). The inset shows the temperature variation of the 
normalized amplitudes for the main frequencies for sample 
B (S = 32°) (data are offset for clarity). Solid lines are fits 
to extract the effective mass (see text) 



would become two distinct ellipses. They fold onto 
each other in the actual Brillouin zone. The in- 
ner electron sheet has mostly dxz (dyz) character 
whereas the outer one has mostly dxy character. 
The inner sheet has significantly less c-axis warp- 
ing than the outer sheet. Close to the zone cen- 
ter are three hole-like sheets. There are two tubes 
running along the c-axis, both sheets have strong 
dxz (dyz) character although the outer one has al- 
most equal c?j,2_j^2 character. These tubes are cut 
through by the dz^ band which forms a squashed 
spherical shape. Hence, we are left with two tubu- 



lar sheets with strong warping near Z and a third 
three dimensional closed cylinder pocket. The po- 
sition of this (1^2 band is strongly dependent on the 
phosphorous atom position (Zp) and also on the 
doping (which wc simulated using the virtual crys- 
tal approximation). 

4-1- Effect of spin- orbit interaction 

A significant effect of the spin-orbit interaction 
is that, as shown in Fig. [Til it splits the two elec- 
tron Fermi surface sheets into two, with a gap in 
/c-space equal to ~1.2% of 27r/a. The bands are 
non-longer degenerate along the X-M (R-A) line. 
The gap is 50meV close to the Fermi level along 
X-M and 40meV along R-A. Another effect is that 
it also increases the gap at the points where the (1^2 
band crosses the hole tube sheets, and decreases the 
warping of these Fermi surface sheets. The changes 
to the Fermi surface are summarized in Fig. 1111 

4-2. Magnetic breakdown 

Given that the gaps are quite small there is a 
question as to whether the electrons may tunnel 
between these Fermi surface sheets when subjected 
to the strong magnetic fields used in our study. 
This effect, which is known as magnetic breakdown, 
has been studied extensively and is found to occur 
when the field is in excess of a 'breakdown field' 
Hbd- As this critical field is crossed the electrons 
would begin to change their orbits so that they fol- 
low the trajectory of the ellipse which would exist 
if there were no gap. There is a smooth crossover 
at Hbd', for H <^ Hbd orbits on the separate in- 
ner and outer electron sheets would be visible and 
for H Hbd there would be just one orbit on the 
joined up ellipse (or two if c-axis warping is taken 
into account). In the intermediate region both fre- 
quencies corresponding to both types of orbit would 
be observed. In general frequencies corresponding 
to Fa + \niFg — Fa), where n=l,2 or 3 might be 
expected [33| . An estimate of the breakdown field 



can be made using the following formula [34 1 



f^oHBD - ^ ^ 



where kg is the fc-space gap between the two orbits 
(measured at the closest point) and and b^^ are 
the fc-space radii of curvatures of the orbits on each 
side of the gap. Using the parameters from our cal- 
culations we estimate iiqHbd — SOT. This would 



increase by a factor ^--^3 to if the real experimen- 
tally determined gaps are used (see later). As Fg ~ 
Fa, + [Fp, - and F, ~ F^, + [Fp, - F„J/2 

these two frequencies could arise from magnetic 
breakdown, although our estimate of Hb d is a fac- 
tor 2 — 6 higher than our maximum field. To rule 
this out completely measurements on higher quality 
samples are required, so that we see if these orbits 
can be observed at much lower field. 

We note here that although the Fermi surfaces 
will be Zeeman split by the strong applied field, the 
measured dHvA frequencies actually correspond to 
to the Fermi surface area Ak{H) — HdAk{H) / dH 
(i.e., the area at a particular field extrapolated to 
zero field with the local slope dAk/dH )|35!]. So 
only non-linear splitting (which we do not expect 
for LaFePO) gives rise to extra dHvA frequencies 
in high field. 

4-3. Comparison to data 

To compare with the experimental results we cal- 
culate the extremal cross-sectional areas of each 
Fermi surface sheet and plot these as a function 
of magnetic field angle as the field is rotated from 
parallel to c to parallel to a. The results shown 
in Fig. 1121 ^^id numerical values are given in Table 
[21 Note that rotations from c towards (110) (i.e., 
at 45° to the a-axis) yield results which are almost 
indistinguishable on the scale of Fig. [T2] (the dif- 
ferences are only clearly apparent as the field ap- 
proaches the a — 6 plane). For each tubular sheet of 
Fermi surface there are two branches correspond- 
ing to the maximal and minimal frequencies. The 
difference in minimum value and angle dependence 
reflects the degree of warping of the tube. 

The calculation for the outer electron sheet (band 
4) is in good agreement with the data for the /3 fre- 
quencies, both in absolute size and warping. The 
inner electron sheet (band 5) is mostly likely associ- 
ated with the a frequencies. The angle dependence 
of band 5 and its warping is in good agreement with 
the data but the absolute value is ~ 30% too large. 
Another reason for believing this assignment is that 
the scattering rate (and absolute signal amplitude) 
for the a frequencies is similar to the (3 frequen- 
cies, which is understandable as both arise from the 
same section of Fermi surface before the spin-orbit 
interaction split them apart. All the a and (3 orbits 
have very similar effective masses which is also un- 
derstandable if they originate from essentially the 
same Fermi surface. The calculations are brought 
into much better agreement with the data if band 4 




Figure 9: The calculated Fermi surface of LaFePO (without 
spin-orbit corrections). 

and 5 are shifted by -|-30meV and -1-85 meV respec- 
tively (see Fig.[T2lJb)). Note that this has the effect 
of increasing the splitting between the two bands 
by a further 32 meV (recall the spin-orbit induced 
splitting was ~ 50meV). Hence, with these shifted 
bands there is now quite sizeable gap in fc-space 
between the two Fermi surface sheets (see Fig.fTTl) 

Without the strongly c-axis dispersive d^'i band 
there would be just two hole-like cylinders at the 
zone centre with very weak c-axis warping. This 
would give rise to two very close frequencies, with 
1/cos 9 angle dependence for each sheet. Indeed, 
this is exactly what is found for the minimum fre- 
quencies of bands 2 and 3. The presence of the 
(i^2 band causes significant warping of each of these 
sheets near their maximum cross-sections with a 
'bulge' around the Z point, which is roughly spher- 
ical locally (see Fig. [TTb). Consequently the angle 
dependence of the maximal frequencies for these 
orbits is very weak, and the maxima have much 
higher frequencies than the minima. The small 
closed three-dimensional cylindrical pocket (band 
1), has a roughly 1/ cos 9 dependence close to 6* = 
and l/sm9 near to 9 ^ 90°. 

With the a and (3 orbits identified, it is natural 
then to identify the weaker 7, 5, and e branches 
with the hole orbits. Given that we do not ob- 
serve experimentally any orbits close to 9 — 90°, 
and the remaining orbits vary like 1 / cos 9 up to 
at least 45°, it seems that the three dimensional 
pocket (band 1) has not been observed. The lowest 
frequency we observe 7 can then be identified with 
the minimum of the inner hole tube. If we shift 
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Figure 11: Fermi surface slices including the spin-orbit interaction (compare to Fig. llOj l (a-d) are (001) slices showing (a,b) the 
hole sheets at r,Z and (c,d) the electron sheets at M,A, (e) shows a (110) slice. The right panels show the same cuts but the 
with bands shifted to best fit the dHvA data (as described in the text) 



all the hole bands by -53 meV we find almost per- 
fect agreement between band 2 and 7, and then the 
minimum of band 3 coincides with d without any 
further adjustment. With this shift the maximum 
of band 2 overlaps with the a orbits so may be dif- 
ficult to observe experimentally. The maximum of 
band 3 coincides with the e orbit although it does 
not have the correct curvature. 

Given that we see no evidence for band 1, and 
that the observed frequencies all have a 1/ cos 6 an- 
gle dependence, so that there is no evidence for 
the 'bulge' close to the Z point (see Fig. [TTfe). it 
is likely that in reality the Fe (1^2 does not actu- 
ally cross the Fermi level As we noted above 



this band is the most sensitive to P position and 
doping. So it could be that 7 corresponds to both 
frequencies from the now weakly warped inner hole 
pocket and i5 (or possibly) e corresponds to both 
frequencies from the now weakly warped outer hole 
pocket. The problem with this is that it leaves e (or 
6) unidentified. Of course, this difficulty would be 
removed if S was really caused by torque interaction 
or magnetic breakdown. Further experiments, with 
lower scattering rates for the hole orbits, as well 
as following how the frequencies move with doping 
and pressure will hopefully resolve these issues. 

From valence considerations we expect LaFePO 
to have exactly equal number of electrons and holes. 
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Table 2: Calculated dHvA frequencies and band masses nii,. Results are shown with and without the rigid band shifts described 
in the text. The results of the fully relativistic band structure calculation reported in Ref. [12|| are also included for comparison 
(denoted S). 
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0.741 
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1.100 


0.741 
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0.547 
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0.992 
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1.1 


0.689 
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1.8 
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0.771 


3Z-h 


3.026 


1.721 


1.170 
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0.868 
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0.919 
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0.7 


0.490 
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0.9 
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0.657 
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0.376 
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1.048 


0.836 


0.9 


0.8 
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Figure 12: Calculated dHvA frequencies versus field angle. The experimental data are shown by black dots, a) shows the 
calculation (including spin-orbit interaction) without any band shifting , and b) with the bands shifted as described in the text. 



In Table [3] we show calculations of the number of 
electron/holes in each of Fermi surface sheets for 
both shifted and unshifted bands. The very slight 



departure of the unshifted bands total from zero 
(0.002 electrons per unit cell) reflects the precision 
of the numerical integration and Fermi surface de- 
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Figure 10: Slices through the Fermi surface calculated with- 
out spin-orbit interaction. The slice are (a) along (001) at 
the centre of the zone, (b) along (001) at the top of the zone 
and (c) along (110) through the centres of both electron and 
hole sheet. The symmetry points are indicated. Each k 
point is coloured according to the maximum band character 
at that point (see key at bottom of figure). 

termination. For the shifted bands the imbalance 
is 0.067 electrons per unit cell. In principle, this 
could suggest that there either remains some sheet 
of Fermi surface which is undetected, or our assign- 



Table 3: Calculated number of electrons (n) or holes in the 
first Brillouin zone of LaFePO for each Fermi surface sheet, 
according to the band-structure calculations (with spin-orbit 
corrections) with (ns) without (njvs) shifting the bands. 
The negative sign on the electron count indicates holes. Note 
this corresponds to two formula units. The shifts (AE) are 
in units of meV. 



Band 




AE 


n-NS 


1 


-0.016 


-53 


-0.001 


2 


-0.098 


-53 


-0.062 


3 


-0.180 


-53 


-0.105 


4 


+0.178 


+30 


+0.162 


5 


+0.118 


+85 


+0.074 


total 


+0.002 




+0.069 



ment of the orbits is incorrect. A possible candi- 
date for the former explanation is the presence of 
the 3D hole pocket which we have not observed - 
however given the likely small size of this it is un- 
likely to make up this difference. Alternatively, if 
we assume that the inner hole pocket is 7 and the 
outer one is e (and attribute S to the torque in- 
teraction) the imbalance would be approximately 
halved, 0.036 electrons per unit cell. This small re- 
maining imbalance could easily be explained by a 
slight oxygen deficiency in our samples. An imbal- 
ance of 0.036 electrons per unit cell corresponds to 
an oxygen deficiency of just 0.9% which is below 
our x-ray detection limit. 

The band structure and dHvA frequencies for 
LaFePO were also calculated by Sugawara et al. 
[l3 | using an in-house fuUy-relativistic code. The 
Fermi surface again is composed for quasi-two di- 
mensional electron and hole sheets however in their 
results the d^^i does not cross the Fermi level. In 
this calculation the sizes of the orbits are all much 
smaller than those found in our work (as well as 
those found experimentally, see TablcJT]). For exam- 
ple, the frequencies for all of the electron orbits are 
less than 1 kT compared to 1.9-2.4 kT found here 
(see Tabled). 

With the various orbits provisionally assigned, we 
can now compare the effective masses. The theo- 
retical values for both shifted and unshifted bands 
are shown in Table [H For most orbits the band 
masses are close to 1. The larger values for some 
of the hole orbits are caused by the top of the (1^2 
passing close to the Fermi level. Comparing with 
the measured effective masses we can see that the 
enhancement factor is approximately 2 for all the 
observed orbits. 



By approximating the Fermi surface sheets as un- 
warped cyhnders, we can estimate the correspon- 
dence between our measured effective masses and 
the electronic specific heat coefficient (je)- In 
the 2D limit, the contribution to 7^ from each 
two-dimensional sheet is given by ^^a^NAkg. 
Four quasi- two-dimensional sheets with m*=2.0 uie 
would therefore correspond to 6 mJ/molK^ (note 
that here a mol refers to a formula unit not the unit 
cell). Experimentally, a value of 'Je = 7.0 ± 0.2 
mJ/mol was estimated for our crystals [l6| 
which is close to our dHvA estimate. 

One remaining question about the bandstruc- 
ture is the influence of the magnetic interactions. 
The experimentally determined positions of the La 
and P atoms are not in their lowest energy posi- 
tions according to the DFT calculations. The cal- 
culated optimal position of the P atom is around 
0.05 A from its experimental position. This error 
is reduced considerably if the ground state is as- 
sumed to be antiferromagnetic '37'| - which it is not 
experimentally found to be. This problem may orig- 
inate from the mean-field nature of the DFT calcu- 
lations which do not take into account the fluctua- 
tion effects which can suppress the magnetism f37| . 
Clearly there are some internal consistency prob- 
lems with the calculations and some authors have 
chosen to base their calculations on the optimized 
atomic positions rather than the experimental ones 
[3^. For the case of LaFePO we find that this ac- 
tually makes the correspondence with the experi- 
mental data worse. Much larger shifts are needed 
to bring the predicted frequencies into agreement 
with experiment. 

We finish by comparing the dHvA results with 
those of a recent angle-resolved photoemission spec- 
troscopy study of the same compound by Lu et 
al. [al]. The ARPES study could follow the 
band dispersions to ~0.5eV below the Fermi level. 
They found that their band-dispersions were in 
agreement with the calculated band-structure when 
renormalized by a factor 2 and shifted by ~ 
lOOmeV. This renormalisation is the same as we 
found for quasiparticles at the Fermi level, imply- 
ing that the renormalisation occurs over the whole 
band, rather than being localized to energies close 
to the Fermi level. Renormalisation over a similar 
wide energy range was also found in Sr2Ru04 [s^. 
The ARPES intensity close to the Fermi level gives 
evidence for small hole pocket (s) at the zone centre 
and electron pocket (s) at the zone corner in agree- 
ment with our results and the band-structure. The 



ARPES results also show evidence for another large 
hole pocket cantered on F. Including this pocket in 
the sum of electron and hole Fermi surface volumes 
leads to a very large electron imbalance of around 
one electron per unit cell, suggesting a significant 
doping of the surface layer. This is possibly caused 
by the polar nature of the cleaved surface [11] . 

5. Conclusions 

Measurements of the de Haas-van Alphen effect 
have been shown to provide a very precise probe of 
the Fermi surface properties of LaFePO. The Fermi 
surface topology deduced from the measurements is 
in good overall agreement with density functional 
theory calculations, consisting of two hole sheets 
near the zone centre and two electron sheets can- 
tered on the zone corner. The data give an accu- 
rate picture of the warping of the electron sheets 
although there remains some uncertainty about the 
correct identification for the hole sheets. We hope 
these issues will be resolved in the near future as 
larger, higher purity samples become available. 
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